On the spectrum behavior of vibrated granular matter 
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Abstract 

A laser facility based on a linear image sensor with a sampling period of 100/is allows to investigate the dissipative dynamics 
of a vibrated granular matter under gravity. The laser reveals the vertical movement of an individual Zirconia-Ytria stabilized 
2mm ball at the surface of a weakly excited 3D granular matter bed. The stochastic realizations are measured from the top 
of the container. Then, power spectra measurements reveal the different cooperative dynamics of the fiuidized gap. We also 
carried out measurements for one steel ball and many balls in ID and 3D systems. We fit the measured different regimes with 
generalized Langevin pictures. We introduce a fractional temporal operator to characterize the ensemble of dissipative particles 
which cannot be represented by a single Langevin particle in a complex fluid. 

PACS: 05.40.-a,47.50.+d,81.05.Rm,83.70.Fm. 



INTRODUCTION 

Granular media (GM) are ubiquitous in nature and 
exhibit a wealth of intriguing physical properties In 
this letter we study the correlations and the statistical 
properties of vibrated GM under gravity. Since energy is 
constantly being added to the system, a non-equilibrium 
steady state is reached |2J, kg • We show that in the weakly 
excited regime the dynamical behavior of the fiuidized 
particles cannot be described as simple Brownian parti- 
cles, this fact leads us to the conclusions that the coop- 
erative dissipative dynamics of the GM particles can be 
described in terms of generalized Langevin particles. 

Hayakawa and Hong Q introduced the thermody- 
namics approach for weakly excited GM. They stud- 
ied vibrated GM under gravity, and mapped the non- 
equilibrium system with a spinless Fermion-like theory. 
The important point is that by focusing on the configu- 
rational properties of an excluded volume (hole) theory, 
the GM non-equilibrium steady state can be understood 
in terms of a configurational (maximum principle) Fermi- 
like distribution. 

Many interesting questions concerning the complex 
system of GM are still open, i.e., the stochastic motion 
z(t) of the fiuidized particles is not entirely known 
For example it is important to test whether the real- 
ization spectrum, S z (f), behaves as a Brownian motion 
(l// 2 ), or it has a more complex behavior related to a co- 
operative dynamics. Therefore an analysis of the stochas- 
tic realization z(t) of these macroscopic Fermi-like parti- 
cles (mFp) should be made. In this letter we report some 
results concerning the analysis of these mFp, in particular 
we have measured these realizations z(t) and from their 
Fourier (FFT) analysis we conclude that S z (f) ~ f~ v 
with v = z/(/ e ,r), (where V = Auj 2 /g, A is the am- 



plitude and u> = 2irf e the frequency of the bed oscilla- 
tions) , showing an agreement with our phenomenological 
(Langevin) memory-like stochastic differential equation 
(SDE) approach. 

A sinusoidal vibration is driven by a vibration plate 
on the GM bed. The vibration apparatus is set up by an 
electromagnetic shaker [TIRAVIB 5212], which allows for 
feedback through a piezoelectric accelerometer (PCB)@. 
This allow us the control of the frequency f e and the 
acceleration in the range of 10Hz to 7000Hz, and 2g to 
40<7 respectively, where g is the acceleration of gravity. 
The control loop is completed by an Oscillator Lab- works 
SC121 and a TIRA 19/z amplifier of lkw. 

The 13 layer GM was set up with Z r C>2 — Y2O3 balls 
of diameter d = 1.99mm into a glass container of 30mm 
of diameter with steel bottom, see Fig. 1(a). Also exper- 
iments with steel balls of d — 10.7mm in different glass 
cylinders (11.2 and 11.4 mm) with steel bottom were car- 
ried out. 

In all cases, the stochastic realizations z{t) in the direc- 
tion of the acceleration of gravity were measured from the 
top of cylindrical containers, during the vibration under 
gravity and at a given f e . The experiments were carried 
out in a chamber at latm of air with 5.8 ± 0.2gr/m 3 of 
water vapor. The absolute humidity was controlled by 
using a peltier condenser and a control loop through a 
thermo- hygrometer. In the case of the experiment with 
Z r Oi — Y2O3 balls, the humidity is of major relevance in 
order to control the contact forces particle-particle and 
particle- wall 0, Under this humidity controlled condi- 
tions, it was not observed surface convection, convection 
rolls in the GM bed, nor rotational movement of the bed 
with respect to the container, which is typical for a con- 
tent of water vapor > 10gr/m 3 . 

The realization z(t) of one particle was followed in 
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FIG. 1: (a) Geometry of the distance measurement [realiza- 
tions of the amplitude z(t)] corresponding to 13 layer GM 
with d — 1.99mm Z r Oi — YiOz balls into a glass container 
of 30mm of diameter, (b) A typical realization z(t) for one 
steel ball (d — 10.7mm) in a ID column, when it is perform- 
ing a quasi non erratic movement under gravity is shown. In 
this case f e = 250Hz and T = 20, thus corresponding to the 
relative length A/d ~ 0.0074 and dimensionless maximum 
velocity V b = Aujj^d ~ 0.39. 



a ID window of 12mm with a laser device by using 

a triangulation method, see Fig. 1(a). In Fig. 1(b) we 
show a typical realization z(t) when one steel ball is 
vibrated at f e = 250Hz with T = 20; under such ex- 
perimental conditions (sinusoidal excitation amplitude 
A = Tg/ut 2 ~ 79.6/zm) it is possible to observe the quasi 
non erratic parabolas for the movement of the ball un- 
der gravity. A laser emitter with a spot of 70/im and 
a linear image sensor (CCD-like array) enables a high 
speed measurement with 100/isec sampling. The linear 
image sensing method measures the peak position values 
for the light spots and suppress the perturbation of sec- 
ondary peaks, which makes possible a resolution of 1/im. 
The shaker and the laser displacement sensor were placed 
on vibration-isolated tables to isolate them from the ex- 
ternal vibrations, and the displacement sensor from the 
experiment vibrations. The realization z(t) is a measure 
of the oscillations of the distance (difference) between the 
particle and the sensor, around the surface of the granu- 
lar bed (the fluidized 3D gap). The measured realization 
without excitation reveals a white noise < 10/mi. Then 
due to the experimental set-up our effective resolution is 
no higher than 10/im. The two-second long register of 

z(t) and the FFT analysis with Nyquist frequencies of 



0.5 and 1.25kHz. were taken with a 9354 C Le Croy Os- 
cilloscope of 500MHz. The velocity V{t) = dz/dt of the 
mFp and its square dispersion Oy = (V(t) 2 } — (V(t)) 2 
were calculated from the z(t) registers. In Fig. 2(a) we re- 



port <7y against a z — <J (z{t) 2 ) — (z(t)) for fixed T = 20 
and several f e from 200Hz to 340Hz for one steel ball ex- 
periments (circles and squares), and also from 250Hz to 
400Hz for nine steel balls in a column. The experiments 
with many Z r C>2 — Y2O3 balls (d = 1.99mm) were carried 
out for fixed T = 20 and f e from 70Hz to 600Hz, thus 
the GM bed was driven by a excitation of amplitude A 
from 1mm to 13/im. Complementary, one minute films 
of this fluidized 3D gap were taking by using a video 
camera microscope (CVM) . For the propose of this work 
the CVM was only used to identify qualitatively the dif- 
ferent dynamical regimes of the Z r C>2 — Y2O3 balls, in 
fact the recorded movies were only used to discuss the 
physical meaning of the measured z(t) registers and their 
corresponding realization spectra S z (f). These measure- 
ments at r = 20 are a partial research of a large work on 
progress we carried out with 5 < T < 40. 



THE SPECTRUM OF THE REALIZATIONS 

When the input of energy is large the individual motion 
of the beads looks erratic, in this section we will analyze 
the frequency spectrum of the stochastic realizations z(t) 
of the mFp at the fluidized gap. We remark that in this 
report we will not be concerned in a regime where the 
input of energy is so small that the power spectrum can 
be described by a Feigenbaun scenario |3j ■ 

Consider a given stochastic process (SP) u(t). Then 
under quite general conditions if the characteristic func- 
tion of the SP u(t) (Fourier transform of the probability 
distribution G u {k 1 1) = [P(u, t)]) obeys the asymptotic 
scaling 



G u (k, t), 



(1) 



which means that the SP u(t) fulfills in distribution the 
scaling A~ H u(At) — > u(t), then, the spectrum density 
behaves asymptotically like 



S u (/)oc l/f 2H+1 . 

Now let us apply this result to a Langevin-like particle. 
Consider the case when the velocity V(t) and the position 
z(t) of the particle are described by the SDE 



dV(t) 
dt 



-lV(t) + t(t), and 



dz 
~dt 



V(t), (2) 



when £(t) is a Gaussian colored noise. This SDE cor- 
responds to the alternative treatment of Brownian mo- 
tion, in fact initiated by Langevin [ToL . This colored 
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noise case corresponds to a weak non-Markovian process 
[HI- From © and (JTJ we obtain, in the overdamped case 
(Brownian motion) S z (f) oc 1/f 2 (limit / < / e « 7). 
On the contrary, in the undamped free motion case (ran- 
dom accelerations model) we obtain S z (f) oc l// 4 (limit 
/ > f e 3> 7). Nevertheless, in the case when the cor- 
relation of the noise £(t) in J2J is of the long-range type 
°c * -9 witn < < 1 (perhaps describing a 
complex fluid surrounding our test particle), the SP z(t) 
turns out to be a strong non-Markovian process i.e., 
there is no Fokker-Planck asymptotic regime. Then we 
get for the spectrum in the overdamped case 

S z (f) oc 1/f- 8 , if dz/alt = f(f), with < 6 < 1, (3) 

and in the undamped case 

S z (f) oc l// 5 " 9 , if d 2 z/dt 2 = £(i), with < 6 < 1. 

(4) 

For the particular situation when the ensemble of 
mFp cannot be represented by a single particle in a 
"fluid" (Langevin-like particle) we should leave the regu- 
lar stochastic calculus and introduce a fractional tempo- 
ral differential operator to emulate the complex multiple 
particle collisions and particle motion. Then we model 
the GM system by a fractional Langevin equation of the 
form 

D? TO-^r^y = -7 a V(t)H(t), and J = V(t) 

(5) 

where £(t) is a Gaussian white noise. So, asymptotically 
in the undamped limit, / > f e 3> 7, we get 

fc 1 
Gv(x^,At) — > G v (k,t), withi? = --(1-a), when a G 

(6) 

which means that Sy(f) oc 1/ f 2a . Then using dzj dt = 
V(t) it follows 

= ^Sy(f) oc 1// 2 ( 1+Q ), a G (i, 1), (7) 

which is the desired result to fit our experiments when 
friction and collisions between the GM particles are im- 
portant issues, and cannot be represented by a memory- 
like Brownian model as in Eq.ljJJ), nor by using the 
memory- like random acceleration model as in Eq. Q . 

The general case given in Eq.JSJ) by keeping arbitrary 
the friction coefficient 7 and for any noise £(t) and cor- 
relation (£(f)£(0)) can also be solved in a similar way 

RESULTS AND DISCUSSION 

In Fig. 2(a) we show the velocity squared dispersion a v 
against the amplitude-dispersion a z of the realizations of 




FIG. 2: (a) cr v as function of a z , the circles correspond to 
one steel ball, full-squares to one steel ball with rotation, and 
triangles to nine steel balls in a ID array. Typical log-log 
power spectrum for one steel ball (b-c) and nine steel balls 
(d) in a ID array. 



the mFp corresponding to the ID array of 10.7mm steel 
balls vibrated vertically; for weak amplitude the slope 
Oyl cr z shows a linear behavior. Fig. 2(a) shows three 
experiments, the first corresponds to one steel ball ver- 
^tic^ly jumping in a cylinder of 11.2mm (almost elastic 
case) with S z {f) ~ l// 4 , see Fig. 2(b), and gives the slope 
(9.7±3)m/s 2 ~ g (in accord with the conservation of en- 
ergy: mgcr z = imcry). The second experiment is one 
steel ball jumping in a wider cylinder (11.4mm), in this 

1 ca^e the slope gives (8.2±3)m/s 2 indicating a dissipative 

2 mechanism, here typically S z (f) ~ 1/f 4 ' 6 see Fig. 2(c). 
The wider cylinder permits the rotation of the steel ball, 
probably due to the interaction with the air gap between 
the steel ball and the wall. In the third experiment there 
are nine steel balls in the wider cylinder, in this case the 
many-body collisions came into account giving a smaller 
slope (5.9 ± 3)m/s 2 , here S z {f) - 1/f 3 - 8 see Fig.2(d) 
(we find that for 400Hz> f e > 200Hz, the exponent is 
3.5 < v < 3.8). This last result indicates that in this 
regime it is necessary to introduce a description in terms 
of a non usual collision operator. The estimated slope 

I Gz gi ves a qualitative magnitude from which it is 
possible to observe a clear tendency in the GM motion 
behavior. Spectral density measurements [corresponding 
to the geometry of Fig. 1(a)] associated to the realization 
z(t) of ZrC>2 — Y2O3 balls are shown in Fig. 3 (a-c). In 
Fig. 3(d) the plot of the exponent v against the excitation 
frequency f e for a fixed acceleration T = 20 is shown. For 
f e = 600Hz, the balls are organized in a "lattice" (we 
test this fact by using a VCM at the surface of the GM 
bed) where the weak- vibration movement of each ball is 
limited to a site of this lattice, given therefore a Wiener 
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FIG. 3: Log-log power spectrum for 13 layers of ZrOi — Y2O3 
balls in a glass container as in Fig. 1(a), for T — 20 and differ- 
ent excitation frequency f e , (a) 73Hz, (b) 137Hz, (c) 225Hz. 
(d) The power spectrum exponent v against the excitation 
frequency f e is shown. 



from the many-body bed. This makes necessary the in- 
troduction of a new collision operator which is written 
in terms of the fractional SDE (JSJ. In other words, the 
hopping of a ball is perturbed (mainly) by few-body col- 
lisions. At frequencies f e around 70Hz, we observe a pure 
collision behavior where a ball is under the random ac- 
celeration mechanism transferred from the surrounding 
"fluid", as it is proved because v is between 4 and 5, 
corresponding to the random acceleration dynamics em- 
ulated by Eq.(U). 

In general our SDE picture allows us to calculate also 
the velocity distribution of the GM, work along this line 
will be reported elsewhere. 
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realization z(t). Note that for T = 20 and f e = 600Hz the 
excitation amplitude is A ~ 13.8^m, so we do not expect 
hopping motion from one to another site. It is interesting 
to remark that even when the figures A/d ~ 0.007 and 
Auo/\fgd ~ 0.37 for the Zirconia-Ytria balls are similar 
to the one from the experiment with only one steel ball 
(for r = 20 and f e — 250Hz), in the present case the 
stochastic behavior is quite different as can be corrobo- 
rated from the value of v ~ 2, when it is compared with 
the one form Fig.2(b-c). For f e approximately between 
500Hz and 200Hz, the "hopping" of the fluidized balls 
(jumping to unoccupied hole-states) are more frequently, 
thus leading to exponents v from 2 to 3. This is a gen- 
eralized Brownian motion from low to high correlation, 
see Eq.©. When v = 3 is reached we observe a pure 
"hopping" mechanism and it is related to a highly non- 
Markovian memory in the SDE (limit 9 — * 0), i.e., the 
long hopping of one ball is associated with the long-range 
noise correlation of the "fluid" . Such a fluid should be 
interpreted as the hopping "events" of the whole gran- 
ular matter bed. For exponents v between 3 and 4 (if 
f e w 200Hz to 150Hz), the motion of a single particle 
evolves from a pure hopping to a more complex behavior 
where the hopping-balls are perturbed by the collisions 
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